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ABSTRACT 

We use a sample of z 3 Lyman Break Galaxies (LBGs) to examine close pair 
clustering statistics in comparison to LCDM-based models of structure formation. 
Samples are selected by matching the LEG number density, Ug, and by matching the 
observed LEG 3-D correlation function of LEGs over the two-halo term region. We 
show that UV-luminosity abundance matching cannot reproduce the observed data, 
but if subhalos are chosen to reproduce the observed clustering of LEGs we are able 
to reproduce the observed LEG pair fraction, (A'c), defined as the average number of 
companions per galaxy. This model suggests an over abundance of LEGs by a factor 
of ~ 5 over those observed, suggesting that only 1 in 5 halos above a fixed mass hosts 
a galaxy with LEG-like UV luminosity detectable via LEG selection techniques. This 
overdensity is in agreement with the results of a Millennium 2 analysis and with the 
discrepancies noted by previous authors using different types of simulations. We find 
a total observable close pair fraction of 23 ± 0.6 per cent (17.7 ± 0.5) per cent using 
a prototypical cylinder radius in our overdense fiducial model and 8.3 ± 0.5 per cent 
(5.6 ± 0.2 per cent) in an abundance matched model (impurity corrected). For the 
matched spectroscopic slit analysis, we find Ncs{R) = 4.3 ± 1.55(1.0 ± 0.2) per cent 
and Ncs = 5.1±0.2 (1.68±0.02) per cent, the average number of companions observed 
serendipitously in randomly aligned spectroscopic slits, for fiducial slits (abundance 
matched) , whereas the observed fraction of serendipitous spectroscopic close pairs is 
4.7±1.5 per cent using the full LEG sample and 7.1±2.3 per cent for a subsample with 
higher signal-to-noise ratio. We conduct the same analysis on a sample of dark matter 
haloes from the Millennium 2 simulation and find similar results. From the results and 
an analysis of the observed LEG 2-D correlation functions, we show that the standard 
method of halo assignment fails to reproduce the break, or up turn, in the LEG close 
pair behavior at small scale (< 20 /i^^kpc physical). To reconcile these discrepancies 
we suggest that a plausible fraction of LEGs in close pairs with lower mass (higher 
density) than our sample experience interaction-induced enhanced star formation that 
boosts their luminosity sufhciently to be detected in observational sample but are not 
included in the abundance matched simulation sample. 

Key words: cosmology: theory, large-scale structure of universe — galaxies: forma- 
tion, evolution, high-redshift, interactions, statistics 



1 INTRODUCTION 

One of the fundamental predictions of a ACDM model of 
the universe is the liierarchical growth of structure. How- 
ever, direct observations of galaxy mergers, and, by exten- 
* E-mail: jberricr@uark.edu sion, statistics on galaxy mergers, are difficult to obtain due 
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to the long time-scales of the galaxy-galaxy merger pro- 
cess. Correlations between galaxy characteristics and their 
environment suggest that interactions play a role in set- 
ting galaxy pro perties such as star formation rate, colour, 
morp hology (e .g. Toomre fc Toom"r3l 19721: Larson fc Tinslev 



1978; Dressier 1980: Postman fc Gelleij 119841 : iBarton et al 
2000 ; Barton Gillesp ie ct al. 20oi )" However, observational 
studies of mergers and interactions can be difficult due to 
the low luminosities of tidal features and the difficulties in 
quantifying galaxy morphologies. At high redshifts, z 1, 
these problems are exacerbated by the decreased apparent 
luminosity and resolution of the gala x ies be ing studied. 

Since the studies of iHolmbergj (|l937l ) close pairs of 
galaxies have provided an important tool for the evaluation 
of galaxy merger rates by providing counts of merger 
candidates, and for theories of galaxy formation due to 
the importance of galaxy-galaxy mergers in galaxy evolu- 
tion. Close galaxy pair counts, or counts of morphologically 
disturbed systems, have not only been used to provide candi- 
dates for galaxy mergers, but have been used in attempts to 
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In 



Berrier et a l. (2006) we present a method to anal- 



yse the close pair fraction of galaxies in a simulation en- 
vironment, with the close pair fraction (N^) defined as the 
number of galaxies in close pairs in a volume of space nor- 
malised by the total number of galaxies in the sample. This 
analysis demonstrates the viability of estimating the observ- 
able close pair fraction in simulations using simple criteria 
to a ssign galaxies to dar k matter haloes. 

iBerrier eHo] l| 2006l) argue that the close luminous com- 
panion count per galaxy does not track the distinct dark 
matter halo merger rate. Instead, it tracks the luminous 
galaxy merger rate. While a direct connection between the 
two has often been assumed, there is a mismatch because 
multiple galaxies may occupy the same host dark matter 
halo. The same arguments apply to morphological identifi- 
cations of merger remnants, which also do not directly probe 
the host dark halo merger rate. This still leaves close galaxy 
pairs as a tracer of galaxy evolution and as a proxy of the 
galaxy merger rate. 

At high redshift, the dense environment and smaller 
fraction of galaxy clusters [where the large velocity 
dispersion prevents many satellite-satellite mergers, e.g. 
iBerrier et all (|2009l 'l] mean that close pairs of galaxies are 
likely to indicate actual mergers, though estimates of the 
t ime scales of these mergers may still be rather large (se e 
e.g.. lKitzbichler fc Whi'^l2008l : iBertone fc Conselicell2009ll . 
As a result, observations of this process for galaxies at high 
redshift are highly desirable for constraining the high red- 
shift galaxy-galaxy merger rate. 

The Lyman break galaxies (LBGs) are star forming 
galax ies efficiently ident ified using colour selection criteria 
(e.g.. ISteidel et al.lll996l ') and comprise a larg e fraction of 
all luminous galaxies at high redshift (e.g., 'Reddy ct al.| 
l2005l : lMarchesini et al. I b007 ). To date, a few thousand LBG 
spectra and tens of thousands of photometric candidates 



have been obtained, making LBGs a useful, well-studied 
population for hig h redshift galaxy spa tial distribution and 
close pair analysis. IConrov et al.l l|2006f ) used subhalo abun- 
da nce matching techn iques (SHAM) such as those used 
in IBerrier et al.l l|2006l ) and here, to calculate angular cor- 
relation functions for LBGs at high redshifts, z = 3, 4. 
This work suggests that abundance matching techniques 
may be used to sample LBG populations and statistics in 
simulations. SHAM has been tested in a variety of situa- 
tions at both low and high redshift and has been shown 
to be a reasonable tool for matching galaxies to popula- 
tions of dark mat ter halos and generating halo mass - stel- 
lar mass relations (jConrov et al. 20061: Vale fc Qstrikerl2006 ; 



Berrier et al.' zOOd; 'Stewart et al.l 



Guo et alJ i 2010; Mostor ct al. 201C|). Because this technique 
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has been suggested, and indeed used, as a probe of LBG 
clustering statistics, it will provide our starting point in 
this analysis. We also explore matching dark matter halo 
and subhalo correlation functions to the observed clustering 
of z '-^ 3 LBGs. T his technique is similar to the work of 
IConrov et al. I (|2008l ) on z 2 star forming galaxies. 

In this paper, we use a numerical Af— body simulation 
with an an alytically generated substructure, adopting the 
approach of lBerrier^ et al.l l|2006l ). to compare close compan- 
ion counts directly to the observed com panion count for our 
sample of z ~ 3 LBGs from the sur vey oflSteidel et al.|(l2003l . 
hereafter S03) and the survey of ICooke et al.l 1 2005 . here- 
after C05). Our pu r pose is to test the s imple and popula,r 
(|Conrov et al.ll2008l : Istewart et aD '2009': 'Si mha et aLll2012l ) 
theory that galaxies live in subhalos and that UV luminosity 
correlates monotonically with halo mass/maximum circular 
velocity at the time of accretion. 

The structure of this paper is as follows. In § [21 we 
outline our methods, discuss our observational sample in 
§ 12.11 our simulations in § 12.21 the models used for the as- 
signment of galaxies to haloes in § 12.31 The definitions of 
the close companion fraction, the photometric companion 
fraction, and the sample impurity and number density are 
covered in § 12.41 - § 12.81 respectively. We present our pre- 
dictions for the companion fraction, Nc, in § O We begin 
with an examination of Nc from 2 = — 3 with an emphasis 
on a comparison between our simulations and the observa- 
tional values at 2 = 3 in § 13.11 The angular photometric 
close companion count is the topic of § 13.21 Comparisons 
with previously existing close companion counts are made 
in § 13.31 We return to the number density issue in § 13.41 
Finally, we discuss the implications of our results in § 13.61 
We conclude with a summary in § 2] 

In this work we assume a flat universe with a standard 
cosmology of flm = 1 — i^A = 0.3, h — 1.0, and as — 0.9. 



2 METHODS 

P air count statistics a re generated using the same technique 
as IBerrier et all l|2006l '). A ACDM A''-body simulation is used 
to identify the large-scale structure and properties of the 
host dark matter halo (details in ^ 12.211 . The analytic sub- 
structure model of IZentner et al.l (|2005l ) is used to generate 
four sets of satellite galaxies within these host haloes for 
our analysis. Using the analytic models with no inherent 
resolution limits to model substructure allows us to over- 
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come the issue of numerical over-merlins in the dense envi- 
ronments (e.g.. [Klvpin et al., 1999). This method has been 
demonstrated to accurately model the two-point clust ering 
statistics of haloes and subhaloes (jZentner et al.l 20051) and 
used to produce viable close pair statistics I Berrier et al.l 
bOOel ) from z = to z = 1. 

We use a simple method to assign galaxies to dark mat- 
ter haloes in our simulation volume (ij |2.3[) and address pos- 
sible effects of this assignment in § 13.61 We conduct mock 
observations on the "galaxy" catalogs in an identical man- 
ner as those used in observational studies to calculate the 
average number of close companions, Nc, or the close pair 
fraction statistic in our simulation box (§ 12. 4p . This can 
be done to mimic the exact specifications of observations 
in the real Universe. The analytic subhalo model allows us 
to examine the variance in close-companion counts associ- 
ated with the realisation-to-realisation scatter. In this way 
we may examine different sets of substructure populations 
while retaining the large-scale structure in our simulation 
allowing us to test for the importance of cosmic variance 
and chance projections. 

In this work we focus on examining the close pair frac- 
tion of potential LBG haloes at 2 = 3 in a simulation box 
and make direct comparisons to sets of observational data. 
In order to more accurately test the expectations of detect- 
ing serendipitous close pairs in conventional multi-object 
spectroscopic surveys, we calculate both a standard Nc, by 
using a cylindrical geometry, and using a mock slit geome- 
try, Ncs, that mimics typical spectroscopic observations and 
those of our survey (ii l2.ip . The mock spectroscopic slits are 
rotated through several possible orientations in the simu- 
lation to calculate the possible variations in observed pair 
fraction caused by the random alignment of the spectro- 
scopic slitlets and the orientation of the galaxy pairs on the 
sky. Our sample of potential LBGs are identified in the simu- 
lation by matching the two point correlation function of ob- 
jects with a given minimum infalling velocity to the observed 
correlation functions. Using lines of sight through the entire 
length of the simulation box, we are able to approximate the 
projected close pair count of LBGs over a defined redshift 
path. Finally, we use multiple randomly aligned copies of the 
simulation box to explore the full line of sight depth of the 
observed sample as a means to test our simulation results 
against the full redshift range of the observations. 

2.1 Observations 

We design certain aspects of the simulation analysis for a 
direct comparison to the imaging and spectroscopic z ~ 3 
LBG surveys of COS and S03. The survey of COS con- 
sists of deep u'BVRI imaging of nine separate fields over 
~465 square arcmi n using the Low Resolution Imagin 
Spectrometer fLRIS jOke et al.lll995l : rMcCarthv et al.lll99. 
on 10-metre Keck I telescope and the Carnegie Observa- 
tories Spectroscop ic Multislit and Imaging Camera (COS- 
MIC; iKells et al.l [1998) on the 5-metre Hale telescope at 
the Palomar Observatory. Approximately 800 photometric 
LBG candidates were selected in a conventional manner that 
uses their u'BVRI colours. The sample contains 211 colour- 
selected, spectroscopically confirmed 2: ~ 3 LBGs with 
m_R < 25.5 and a redshift distribution of (z) = 3.02, la = 
0.3. The nine fields of the survey minimize the effects of 



cosmic variance. Detailed information regarding the colour- 
selection technique and survey specifics can be found in COS. 
The survey of S03 consists of the publicly available photo- 
metric catalog of ~ 2500 2^3 LBGs and the the spatial 
correlation results using a spectroscopic subsample of ~ 800 
LBGs. 

Although the sensitivity limits of 8-metre class tele- 
scopes enable photometric detection of z ~ 3 LBGs to 
m_R < 27, spectroscopic confirmation is limited to those with 
m_R < 25.5 using reasonable integration times. The spatial 
distribution, or clustering, of the spectroscopic sample has 
been used to infer the average mass of LBGs in the context 
of AC DM cosmology (|Adelberger et al.l l2005l : ICooke et al] 
l2006bl . hereafter COB) and is determined from the m_R < 25.5 
subsample. For comparison to our simulation, we only con- 
sider LBGs that have mn < 25.5 in order to compile a sam- 
ple with (1) accurate photometry (< 0.2 magnitude uncer- 
tainties), (2) follow-up spectroscopic confirmation, and (3) 
a measured spatial correlation function. 

The COS survey is a conventional multi-object spectro- 
scopic (MOS) survey originally designed to obtain a large 
number of z ~ 3 LBG spectra to cross-correlate with quasar 
absorption line systems. Although it is unclear whether the 
presence of quasars in these fields produces a clustering bias 
for LBGs near the same redshift range, the background 
quasars for six of the nine fields surveyed are at a much 
higher redshifts than the (3.0), Isigma = 0.3 LBGs probed 
(see COS), thus eliminating any potential clustering bias. 
Any clustering bias for the remaining three fields is likely 
small because the LBG correlation values for the nine fields 
in our survey agree, within the u ncertainties, to the results 
of lAdelberger et all (|200l 120051 ) on the 17 field survey of 
S03. Nevertheless, we generate our simulation sample based 
on the values of S03 to help alleviate any potential bias. Fi- 
nally, we note that two of the serendipitous spectroscopic 
close pairs in our survey are found in the three fields poten- 
tially biased by the targeted quasars but exist at much differ- 
ent redshifts as compared to the quasars {5^ corresponding 
to > 200 h"^ Mpc, comoving) as not to be biased. 

Conventional MOS surveys of LBGs target single LBGs, 
not LBG pairs, and are designed to typically have the same 
orientation for the multiple slitlets located on each slitmask. 
As such, the slitlets have orientations that are random with 
respect to the orientation of LBG pairs on the sky. As a 
result, the fraction of serendipitous LBG pairs that fall into 
the the MOS slitlets enables an accurate sampling of the 
underlying close pair fraction. An illustration of this con- 
cept for one of the many slitlets on a multiobject slitmask 
is shown in Figure [1] Although LBGs cluster, the relative 
low surface density of z ~ 3 LBGs resu lts in very few pairs 
falling serendipitously into the slitlets. ICooke et all (|201Gl . 
hereafter CIO) identify 10 LBGs in 5 serendipitous spectro- 
scopic close pairs (< 20h~^ kpc, physical). The serendipitous 
close pairs provide spectroscopically identified interacting 
events to compliment photometric close pairs and morpho- 
logical classifications which have previously been the only 
means to identify high redshift interactions. Finally, because 
the instruments, method, and analysis of our survey are vir- 
tually identical to most other conventionally acquired sur- 
veys, and specifically to that of S03, it is valid to compare 
the overall results from this work. 

Typical z ~ 3 LBG spectra have a signal-to-noise ratio 
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Figure 1. Illustration of conventional s pectr oscopic slit geometry 
(data from the survey of lCooke et alj ||2005| )) and mocked- up in 
our simulation analysis. In both panels, the geometry of the spec- 
troscopic slitlets is shown by the rectangles, the colour-selected 
2 ~ 3 LBGs are marked using small circles, and a radius 30h~^ 
kpc at z = 3 is denoted by a large (dashed) circle centred on the 
targeted LBG. Although the actual slitlet dimensions vary for 
each target, those illustrated here have the average dimensions of 
our survey (see text). The direction of the spectroscopic dispersion 
precludes acquisition of objects to the left and right of the slitlets 
as depicted in the two panels. For example in the left panel, the 
close companion to the immediate right of the targeted galaxy, 
as well as the more distant companion to the upper right in this 
highly clustered case, cannot receive spectroscopy and must await 
future observations which are not always possible. Occasionally, 
the bulk of the flux of an LBG pair will fall serendipitously into 
a slitlet as shown in the right panel. Each slitlet in the survey is 
among ~ 30 — 40 similarly aligned slitlets acquired per telescope 
pointing that are oriented to minimize atmospheric dispersion at 
the time of the observations and not designed to align with the 
orientation of close pairs. As a result, any LBG pairs that fall 
into the slitlets randomly sample the true underlying close pair 
fraction. 



(S/N) of only a few, but in practice the strong UV emission 
and absorption features and continuum profiles provide a 
means for reliable identification. Nevertheless, cautious of 
the inherent low S/N, we assign a confidence qualifier to the 
spectroscopic identifications. For our pair analysis, we test 
two samples from the observations: The full sample of 211 
LBGs and a sample of 140 LBGs with the highest S/N which 

we term the highest confidence sample. 

The color-sele ction technique ((e.g. ISteidel et al.| [20o3 : 
ICooke et al.) |2005| ) is highly efhcient in targeting 2: ~ 3 
LBGs and removing background and foreground sources. 
The observed 2-D color-selected close pair fractions were 
estimated after a correction for chance alignments by gener- 
ating random catalogs matched to the density, dimensions, 
and photome t ric se lection functions specific to the C05 and 
ISteidel et all (|2003l) surveyed fields. 

2.2 Simulations 

The simulation used for the large scale structure and host 
haloes was performed using the Adaptive R efinement Tree 
(ART) Af-body code (|Kravtsov et al.lll997l ) for a universe 
with a standard cosmology of Qm = 1 — = 0.3, h — 0.7, 
and as = 0.9. The simulation followed the evolution of 512^ 
particles in a comoving box of 120/i"^Mpc on a side, with a 
particle m ass of nip ~ 1.07 x l O^fe' ^M^. More det a ils ca n 
be found in lAllgood eFal] (|2006l ) and lWechsler et all (|2006l ). 
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Figure 2. The cumulative mass function of host haloes as de- 
rived from our 120 h~^Mpc simulation box at 2 = 3. Error bars 
estimate cosmic variance using jackknife errors from the eight 
octants of the computational volume. 
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Figure 3. The vertical axis shows the cumulative number den- 
sity of galaxies in our simulation catalog as a function of velocity 
using Vin (solid black line) and Vnow (blue dashed line) at 2 = 3 
to identify subhaloes as galaxies . The error bars shown were 
generated by summing in quadrature the jackknife error and the 
realisation-to-realisation scatter and represents errors due to cos- 
mic variance in the simulation. 
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The root computational grid was comprised of 512^ cells 
and was adaptively refined according to the evolving local 
density field to a maximum of 8 levels. The peak spatial 
resolution is /ipoak — 1-8 h~^kpc in comoving units. 

In this simulation the distinct host haloes are identified 
using a variation of the B ound Density Maxima algorithm 
(BDM. lKh^pin et al.|[l999l ). In this method each halo is asso- 
ciated with a density peak. This peak is identified using the 
density field smooth ed with a 24-particle SPH kernel (see 
iKravtsov et al.l [20041 . for details). The halo virial radii and 
mass are calculated for the host halo in the simulation box. 

The halo virial radius, Rvir, is defined as the radius of 
a sphere whose centre is the density peak, with mean den- 
sity Avir(z) times the mean density of the universe. The 
virial overdensity Avir(z) comes from the spherical top- hat 
collapse approximation. In our case this is c o mput ed us- 
ing the fitting function of iBrvan &: Normaj] l| 19981 ). The 
simulation assumes a conventional ACDM cosmology which 
yields A^ir{z = 0) ~ 337 and Avir(2) 178 at z> 1. The 
virial mass is used to characterise the masses of distinct host 
haloes, the haloes whose centres do not lie within the virial 
radius of a larger system. 

Figure [2] shows the host halo mass at 2 = 3 from the 
procedure described above. The Figure illustrates host halo 
mass function, complete to virial masses 
The uncertainties are calculated by a jackknife error tech- 
nique. They are computed by removing one of the eight oc- 
tants of the simulation volume and recalculating the mass 
function. These error bars estimate the uncertainty in host 
halo counts from cosmic variance. 

The substructures originally located in these host haloes 
are removed and replaced by subst r ucture s generated us- 
ing the algorithm of Zcnt ner et al.l l|2005l ). This analytic 
method allows the generation and examination of substruc- 
ture with effectively unlimited resolution. Each host halo in 
this simulation catalog has a randomly generated mass ac- 
creti on history using t he extended P ress-Schechter formal- 
ism (iBond et al.lll99ll: iLacev fc Colelll993 ') with the imple- 
mentation of Somerville fc Kolatj 1 19991 ). 

Once these mass accretion histories and merger trees 
are generated, we track the history of the new subhaloes 
as they evolve. As each subhalo merges into the host it is 
assigned an initial orbital energy and angular momentum. 
Then the routine calculates the orbit of the subhalo inside 
a potential from the host halo between the time of subhalo 
accretion to the epoch of observation. Tidal mass loss and 
dynamical friction are modeled to determine the effects of 
these interactions on the mass of the subhalo. The halo's 
density profile is modeled using the the lNavarro et ah! l|l997l . 
NFW) profile wit h halo con centrations set according to the 
algorithm of Wec hsler et all (2002). Finally, all subhaloes are 
tracked until their maximum circular velocities drop below 
Vuiax = 80 km s~^. Haloes which fall below this threshold 
are removed from the simulation. This step is to avoid excess 
computing time calculating the small, tightly bound orbits 
of objects that are not likely to host a luminous galaxy. We 
refer the interested reader to § 3 of Z05 for the full details 
of this model. 

This process is repeated four times for each host halo to 
determine the effects of variation in the subhalo populations 
with a fixed host halo population. In addition to generating 
substructure catalogs for each separate "realisation" of the 



model we perform three rotations of each simulation vol- 
ume. These rotations provide us with different lines of sight 
through the substructure of the simulation. This provides a 
total of 12 effective realisations for us to gather close pair 
statistics. 

Z05 demonstrates that this method is successful at re- 
producing subhalo count statistics, radial distributions, and 
two-point clustering statistics measured in high-resolution 
A'^-body simulations in the regimes we use here. This model's 
results agree with numerical treatments over 3 orders of 
magnitude, or more, in host halo mass as well as a function 
of redshift. Moreover, this technique has also proved useful 
in generating close galaxy pair cou nts that match obse rva- 
tions in the local universe to z ~ 1 (|Berrier et aLlbOOd l. 

In addition to our primary simulation sample de- 
scribed above, a subsample of dark matter haloes from the 
Millennium-2 simulation (see lBovlan-Kolchin et al ]|2009l , for 
more details) are used to test our methods in a pure A'^-body 
simulation. 



2.3 Assigning Galaxies to Haloes and Subhaloes 

After computing the properties of haloes and subhaloes in 
a ACDM cosmology, the next step is to map galaxies on to 
these objects. We use the maximum circular velocity that 
the subhalo had at the time it was accreted into the host 
halo, Vin, to define the objects in our sample. The choice 
of Vin mimics a case where a galaxy is highly resistant to 
baryonic mass loss when compared to its dark matter halo. 
As such, the luminosity of the galaxy is unchanged by the 
loss of matter due to tidal interactions. This case assumes 
a model in which the luminosity of a galaxy is set in the 
field and does not change after merging into the host sys- 
tem. Thus we assume that there is a monotonic relationship 
between halo circular velocity, Knax, and galaxy luminosity. 
This model does not account for any effects which might 
alter the galaxies intrinsic luminosity or which might inter- 
fere with observations, such as galaxy-galaxy interactions 
triggering enhanced star formation or dust obscuration. In 
effect, this model assumes a perfectly observable universe 
with a strong halo assumption that galaxy properties are 
set by the da rk matter halos they reside in. The results of 
IConrov et al.l (j2006) suggest that this form of SHAM may be 
used to examine the clustering statistics of 2 = 2 — 4 LBGs. 
Recent works have suggested that it is reasonable to assume 
a halo-UV luminosity relation (essentially a hal o-star forma- 
tion r ate relation) at the r edshifts we examine (jSimha et al.l 
l2012l : IConrov et al.ll2008l : IStewart et al.ll2009D . We discuss 
the effects of this method of halo assignment on the results 
in §1311 

In addition to the Vin model, a second toy model is 
tested that uses the Vmax of all haloes at the epoch they are 
observed. We refer to this model as the V„ow model. This 
model describes a physical scenario in which the dark matter 
and luminous baryonic matter are stripped from subhaloes 
proportionally. This is in stark contrast to the Vi„ model 
where the luminous baryons are resistant to mass loss. This 
second model has proved to be inadequate in reproducing 
close pairs of galaxies and features observed in the two point 
correlation function of galaxies locally, but is tested here for 
the purposes of completeness. Although baryons are likely 
to be stripped from a halo, it is unlikely that they will be 
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stripped at the same rate as the dark matter. Again, this 
model does not account for the possibility of enhanced star 
formation due to galaxy interactions. 

Figure [3] shows the cumulative number density of 
"galaxies" identified in our simulations, a function 

of their maximum circular velocities. The black solid line 
shows z = 3 galaxies using V^n as an identifier, while blue 
dashed line uses V„ow to generate the function. Our catalogs 
are complete to a Vmax = 100 km s~^. 

In addition to testing a standard SHAM sample and 
in order to make as direct a comparison as possible with 
observational data, we match the two-point spatial corre- 
lation functions of our model galaxy catalogs to the ob- 
se rved z ~ 3 LBG two-p oint spatial correlation function 
of lAdelberger et all (|2003i ) . The correlation functions are 
shown in Figure [l] We note that the low-resolution spec- 
tra and intrinsic star forming processes of « ^ 3 LBGs make 
it difficult to obtain precise redshifts from the emission and 
absorption features l|Shaplev et al.1 l2003l : lAdelberger et al.l 
[2003 ). We consider the effect of LBG redshift uncertainties 
on the pair fracti ons in the next section. T he spatial corre- 
lation functions of lAdelberger et all l|2003h used here, incor- 
porate LBG angular information and adopt a prescription 
(see Appendix C of that work) that aims to minimize red- 
shifts errors in order to estimate the true three-dimensional 
correlation function. We fit the region from approximately 
the innner separation radius computed by that prescription 
out to higher radii and thus heavily dominated by the two 
halo term region. 

The resulting three dimensional two-point correlation 
function from the simulations is averaged over all four cat- 
alogs, with errors including realisation-to-realisation scat- 
ter and jack knife errors. Both the fit to the observed 
"real space" correlation function and the simulation follow 
a power law of the form: 

C(0 = ir/ro)-\ (1) 

where ro is the spatial correlation length and 7 is the power 
law slope. The analysis of C06 places the correlation length 
at rp = 3.3 ±0.6 using a fixed 7 = 1.6, whereas the results of 
lAdelberger et al.l (|2003l l find a value of n, = 4.0 ± 0.6 with 
7 = 1.57 ± 0.14 for the larger S03 survey dataset. When 
matching the two-point correlation function of the simula- 
tion to the data we find that haloes with Vin ^ 133 km 
s~^ best-fit the parameters of the observations. These haloes 
produce values of ro = 3.93 ± 0.61 and 7 = 1.57 ± 0.05. The 
Vnow model that best matches the observed correlation func- 
tion is found to have Vnow 5S 142 km s~^with ro = 3.99±0.64 
and 7 = 1.54 ± 0.05. 

As discussed above, the spatial distribution, or cluster- 
ing, has been used to infer the a verage mass of LBGs in 
the context of ACDM cosmology. lAdelberger et al.l (120051 1 
and C06 find the mass of the m_R < 25.5 LBG spectroscopic 
sample to be (M) ~ 1O""±°-^M0. The mean mass of our 
sample compares weU at (M) = lO^ '^^/i'^Mo. 

2.4 Defining Close Galaxy Pairs 

Now that we have selected our candidate LBG haloes we 
must determine the best way to "observe" our sample in 
order to make direct comparisons with observational results. 
We examine spectroscopically discernible LBG pairs first. 




r (h~' Mpc) 

Figure 4. The three dimensional two-point correlation func- 
tions. The blue dashed hne, and square points, are the real space 
two point co rrelation function measu red from 2 ~ 3 LBG ob- 
servations of lAdelberger et al. The red dot-dashed line, 
and triangular points, are the real space two poin t correlation 
functio n measured from z ~ 3 LBG observations ICooke et al.l 
ll2006bl) . Both sets of observations are converted from angular 
correlatio n functions to thr e e dim ensions using the approxima- 
tion from lAdelberger et al ] l|2003h . Here the /3 and Ix are the 
Beta function and the incomplete Beta function. The solid black 
line is the spatial two point function generated from the four 
simulations. The uncertainties are generated by a combination 
of jack-knife errors and realisation to realisation scatter. Using 
haloes with Vin 1? 133 km s~^ in our simulations results in the 
best match to the observations. 



We define a spectroscopic close pair three ways. Our 
first criteria include galaxies with a separation of 10 — 
30 /i~^kpc on the sky a relative velocity difference in the 
range of —500 5S Vdiff ^ 500 km s~^. The separation range 
on the sky reflects a conventionally determined distance in 
which close galaxy pairs are considered merger candidates 
and is designed to exclude close pairs that would likely ap- 
pear as a single galaxy in the images. The velocity differ- 
ence, if assumed to be dominated by the peculiar velocities 
of the galaxies, corresponds to the haloes that have a high 
probability of merging. The 10 — 30 /i^^kpc separation is 
well measured by our data. The compact, near point source 
nature of 2: ~ 3 LBGs and the typical FWHM seeing of 
the images allows separation of individual galaxies down to 
~ 6/i~^ kpc and the average length of the slitlets used in 
the spectroscopic observations is ~ 37.7 /i^^kpc at 2 = 3, 
with the bulk of the slitlets probing beyond the conventional 
maximum separation. 

The spectroscopic FWHM resolution of our observa- 
tions is ~ 400 km s~^ but velocity offsets as small as ~ 200 
km s~^ can be measured using multiple cross-correlated 
spectral lines and high-signiflcance Lya features. As a re- 
sult, we test a second criteria that includes pairs with 
0—10 ft^^kpc separations in projection on the sky for ob- 
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jects that may appear indistinguishable in the images but 
show clear indications of two separate spectra with a suffi- 
ciently large velocity difference (Vdiff ~ ±200 to ±500 km 
s-i). 

Finally, we exploit the behaviour of the prominent Lya 
feature in z ~ 3 LBGs which can be observed in absorption, 
emission, or a combination of both. The peak of this feature 
in emission is observed to be redshifted from the sy stemic 
redshift by 450 ± 300 km s~^ l|Adelberger et alj [20031 '). with 
the tail of the distribution extending beyond 1000 km s~^. 
The observed redshifted peaks are attributed to galactic- 
scale outflows driven by stellar and supernovae winds. In 
this picture, the blue wing of the Lya emission feature is 
absorbed by neutral gas moving toward the observer and 
Lya photons traveling away from the observer are shifted 
off-resonance as they scatter off receding portions of the out- 
flow, enabling their escape back toward the observer. 

CIO find that every LBG in the serendipitous spectro- 
scopic close pairs in their sample and every spectroscopic 
LBG with a colour-selected close (< 20 h'^kpc) LBG ex- 
hibits Lya in emission. Because the Lya feature is typically 
detected at high significance, we test a third criteria that 
includes close pairs with no minimum separation on the sky 
and with no Vdiff - The large velocity dispersion Lya peaks 
can help to enable the identification of LBG pairs with lit- 
tle or no actual velocity difference. Random samplings of 
the Lya emission velocity offsets show that nearly all such 
galaxy pairs should be spectroscopically discernible. 

To summarise we test three different criteria to select 
close pairs. Each of these criteria use a maximum outer ra- 
dius of 30 /i~^kpc and a maximum velocity difference of 
±500 km s"\ The remaining parameters for the different 
criteria are: 

• (A) Pairs with minimum separations less than 
10 h-^kpc are always excluded (our fiducial sample for the 
cylinders). This is designated as Nc for cylinders and Ncs{A) 
for spectroscopic slits. 

• (B) Pairs with minimum separations between — 
10 /i~^kpc are included (thus pairs with separations — 
30/i"^ kpc are considered) if there velocity difference is 
Vdiff ^ 200. These results are labeled Nc{B) for cylinders 
and Ncs{B) for spectroscopic slits. This case allows us to 
test an intermediate case between criteria (A) and (C). 

• (C) No minimum separation or velocity difference. All 
pairs with separations < 30h~^ kpc and —500 < Vdiff < 500 
km s~^ are identified. This will be our fiducial sample for the 
spectroscopic slit measurements, Ncs- These are reported as 
Nc{C) for cylinders and Nc for spectroscopic slits. 

The differences in the results of these three criteria are small 
(see §[33]). 

With these parameters we calculate the close pair frac- 
tion of galaxies. This quantity is defined as: 



Here rip is the number density of pairs and Ug is the number 
density of galaxies in the sample volume. Thus Nc reflects 
the fraction of galaxies that have close companions. 



2.5 Spectroscopic Companions 

We can perform an analysis of the underlying close pair frac- 
tion using the serendipitous spectroscopic LBG pairs of CIO 
and by mimicking the observation approach of these data in 
the simulation. In order to best determine the probability 
of observing a serendipitous spectroscopic pair in the sim- 
ulation, we generate slits in two specific ways. Firstly, we 
construct mock slits having the actual lengths and widths 
used in the observations with the objects placed at the ob- 
served locations in the slitlets. We place a randomly assigned 
slit on the candidate LBG haloes in the simulation and then 
rotate the slit through 360 degrees, in steps of 20, to com- 
pute the average Ncs for all pairs "observed" . Within the slit 
geometry we utilise criteria (C) above as our fiducial sam- 
ple. Measurements made in these randomly assigned slits 
are designated A^cs(-R). 

Secondly, we recompute Nca from our simulations using 
a prototypical slit length and width to make a generalized 
"observation". The prototypical slit has a width of 1.37 arc- 
sec, the mean width of the slits used in the observations, 
for a half width of 3.7 h~^kpc. The prototypical slit length 
is 60.0 /i~^kpc, which is 11.13 arcsec at z = 3. All lengths 
are given in proper, physical units. The Ncs measured in the 
spectroscopic slit only incorporates pairs observed in one of 
the rotations of the slit. Thus the Ncs is smaller then the 
total Nc for any of the criteria discussed above. Because 
the candidate LBG haloes are centred in the prototypical 
slits we only need to perform 10 rotations over 180 degrees 
to determine the random slit count. Again, we use criteria 
(C) above as our fiducial sample. Modeling the observations 
in these two ways enables a direct comparison to the CIO 
analysis and provides results that can be applied in a more 
general way to any conventionally acquired survey. 

We may also characterise the dependence of the sam- 
ple on the maximum pair separation used. As observed in 
Figure [5] increasing the maximum radius used for our close 
pair counts increases our resulting Nc- Though the increase 
is comparatively small in the spectroscopic slit, we can see 
in Figure [5] that it is significant for the full cylinder sample. 

2.6 Photometric Companions 

We make one other measurement in this work. We examine 
the apparent angular, or "photometric", close pair fraction, 
Np, and it's impurity. A photometric close pair is defined as 
one in which we have no cut on velocity, and also no pairs 
are allowed within the minimum radius. This is essentially 
a single line of sight through the box designed to mimic 
the companion counts observed in the plane of the sky in 
photometric surveys that utilise simple LBG colour selection 
criteria at z ~ 3. 

Our simulation identifies close pairs within a redshift 
range of 5z ~ 0.18 at z = 3. The criteria of S03 and COS 
select LBG populations with (z) = 3.0, la = 0.3. While our 
simulation samples ~ 25 per cent of the LBGs detectable in 
z ^ 3 surveys, it probes a large enough redshift path to dis- 
tinguish objects that are, and are not, physically clustered. 
This is true in part because clustering effects are negligible 
beyond a radius of ~ 10/i~^Mpc. Thus we can use our anal- 
ysis to provide an estimate of Nc and the sample impurity. 
Our results for this work are found in section 13.21 
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Figure 5. The dependence of the close pair fraction on the 
maximum separation between the central galaxy and its compan- 
ion. All data is calculated at z = 3. The solid black line is the full 
cylinder sample. The blue dashed line is the close pair fraction for 
the spectroscopic slit sample. Uncertainties are calculated from a 
combination of the standard deviations of the four realisations 
and jack knife errors in the individual simulation boxes. 
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Figure 6. The measured sample impurity using a varying outer 
radius. The solid black line shows the impurity fraction for the full 
cylinder sample, while the blue dashed line shows the impurity 
for the randomly oriented spectroscopic slit sample. Uncertainties 
are calculated from a combination of the standard deviations of 
the realisations and jack knife errors in the individual simulation 
boxes. 



2.7 Diagnosing the effects of Interlopers 

Our method characterises the probability of chance projec- 
tions being identified as a companion galaxy. We define the 
sample impurity as the fraction of "observed" close pairs 
in the simulation, using the criteria described above, that 
do not reside inside a mutual dark matter halo and are not 
a physically interacting pair. The total sample impurity is 
given by 



Here rif is the number density of false galaxy pairs in the 
sample volume and ric is the number density of galaxy pairs 
observed in the sample. 

Figure[6]identifies the sample impurity as it evolves with 
the radii of the cylinder used. In this case we see how im- 
purity is effected by the maximum size of the cylinder. This 
is of course a simple relationship. As we increase our max- 
imum radius we have a greater chance of identifying a pair 
of companion galaxies, but also a greater risk of picking up 
a chance projection of two physically unassociated galaxies. 
In this figure the uncertainties are calculated by summing 
in quadrature an error associated with cosmic variance, cal- 
culated by a jack knife error method, and a realisation-to- 
realisation scatter. 

We also examine the effect that extending the maxi- 
mum velocity between galaxies has on both the pair fraction 
and the sample purity. Figure [7] shows this relationship. As 
before, using a larger velocity cut increases both our pair 
fraction and sample impurity. The error bars are calculated 
in the same manner as Figure [G] 



In the observational samples of S03 and C05, the color- 
selection techniques are highly efficient (>90% effective) 
in targeting z ~ 3 LBGs and removing background and 
foreground sources as determined by the spectra. The ob- 
served 2-D color-selected close pair fractions reported below 
were corrected for impurity by generating random catalogs 
matched to the density, dimensions, and photometric selec- 
tion functions specific to the C05 and SOS surveyed fields 
and computing the fraction of random unassociated pairs 
occurring within the projected separations for the criteria 
described above. 



2.8 Galaxy Number Density Problem 

For our z = 3 comparisons we adopted a sample of poten- 
tial LBGs that best describes the data and has Vin ^ 133 
km s~^, that is, the Vmax of a subhalo when it is initially 
accreted into a host, or simply the Vmax of the object if 
the halo is an independent host halo. In addition to match- 
ing the two point correlation function, these models pro- 
duce a comoving number density of galaxies, Ug, that may 
also be compared to the observational sample. The volu- 
metric comoving number densities for the two models are 
ng{Vi„ ^ 133 km s"^) = 0.019 ± 0.003 h^ Mpc"^ and 
ng{V„ow ^ 142 km s-i) = 0.01442 ± 0.00003 h^ Mpc'^ 

The comoving number densities for more massive haloes 
are: ng{Vin > 200 km s"^) = 0.0039 ±0.0006, h^ Mpc"^ and 
ng{Vir^ > 300 kms-1) = 0.00059±0.00010 h^ Mpc"^ . While 
we see that the Vin ^ 200 sample has a better match to 
the number density of ~ 0.004 ± 0.002 h^ Mpc"^ from LEG 
observations (see below for further details) , this sample does 
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Figure 7. The measured sample impurity using various velocity 
cuts. The solid black line indicates the impurity of the sample 
for the full cylinder counts. The dashed blue line is the impu- 
rity recorded for the randomly oriented spectroscopic slit sample. 
Uncertainties are calculated from a combination of the standard 
deviations of the four realisations and jack knife errors in the 
individual simulation boxes. Interestingly this demonstrates that 
large maximum velocities do not significantly increase sample im- 
purity over the range shown here. 



not match the observed correlation function and produces a 
significantly lower value for Nc- 

The matched V„ow sample has a surface density in 
comoving coordinates of 6.14 h^ Mpc~^, 10.2 galaxies ar- 
cminute compared to 8.2 h^ Mpc"'^, 13.6 galaxies ar- 
cminute ~^ for the Vi„ sample. Both of these values ex- 
ceed the ~ 1.7 galaxies arcminute from observation. 
These surface densities are measured from the length of the 
120/i~^Mpc box at z = 3, a redshift range of Sz ^ 0.2 and 
then corrected for the total redshift pathlength observed in 
the surveys and the efficiency of LBG detection as a function 
of redshift from the colour selection (assumed to be 100 per 
cent efficient near 2 = 3, e.g., no lost galaxies as a result of 
bright stars, low S/N regions on the chip, colour detection 
efficiency, etc., as is the case for the observations). 

Previous research has uncovered similar discrepancies 
between the number density of matched massive haloes 
and that observed for z 3 LBGs u s ing different types 



of simulations 



e^ 



Dave et al 



Lacev et al.l 



200C 



lOuchi et all |2004 
201H ). Nevertheless, each 



iNagamine et al.H2005l : 

propose that the excess Ug can be resolved by including other 
types of high redshift objects, such as dust obscured and/or 
low star formation rate galaxies and damped Lyg ab sorption 
systems (DLAs). The model of iLacev et al.l (|201ll ') suggest 
that without dust extinction the number of LBGs would be 
~ 5 times the number observed. This model also recreates 
the properties of the observed sample once dust extinction 
is included. Thus, we report our results using the full (high- 



3 RESULTS 

Using the methods described in Section 12.31 we calculate: 
(1) the close pair fraction observed serendipitously in the 
spectroscopic slits, Ncs, (2) the 10 — 30 h~^kpc total spec- 
troscopic close pair fraction, A'^c, (3) the 10 — 30 h^^kpc 
photometric close pair fraction, Np. 

As a reminder, all three sample criteria include pairs 
with separations < 30 /i~^kpc and a velocity difference of 
±500 km s~^. Ifowever, the differences are that criteria (A) 
excludes all pairs with separations < 10 h'^kpc , criteria (B) 
includes pairs with separations < 10 /i~^kpc if their velocity 
difference is > 200 km s~^, and criteria (C) includes pairs 
with separations of < 10 /i~^kpc with no restriction on the 
velocity difference. 



3.1 Spectroscopic Close Companion Counts at 
z=3 

We first note that the observations of CIO, to which we are 
comparing these results, find a serendipitous close pair frac- 
tion of Ncs = 0.071 ± 0.023 in the spectroscopic slits for 
the highest confidence subsample (140 LBGs) and Ncs ~ 
0.047±0.015 for the full sample of 211 LBGs. In our simula- 
tions we "observe" serendipitous close pairs using the three 
different sets of criteria discussed in § 12.41 

In our randomly selected slit length sample (lengths 
and object positions matched to the survey of C05) we find 
Ncs{R){A) = 0.0220 ±0.0234, Ncs{R){B) = 0.0300±0.0178, 
and Ncs{R)iC) = 0.0430 ± 0.0145. AU three resuhs are con- 
sistent within their uncertainties. 

For the prototypical slit length of 60.0 /i^^kpc, we find 
an "observed" Ncs{A) = 0.0293±0.0013, Ncs{B) = 0.0375± 
0.0016, and Ncs{C) = 0.0506 ± 0.0017. 

As before, the uncertainties are a combination of 
jack knife errors from cosmic variance and realisation-to- 
realisation scatter. These measurements refiect our uncor- 
rected sample and include galaxies that meet our velocity 
criteria and are observed as close pairs due to projection 
effects even though they do not reside in the same parent 
halo. 

Finally, the full observable pair fraction in the cylin- 
der is Nc{A) = 0.228 ± 0.006, Nc{B) = 0.246 ± 0.006, and 
Nc{C) = 0.275 ± 0.007. As expected, criteria (B) and (C) 
produce slightly higher values, but they are not significantly 
different from our fiducial value (A). 

We report the results of criteria (A) for the cylinder and 
criteria (C) for the spectroscopic slits because criteria (A) is 
designed to match the morphological and close pair fraction 
analyses in the literature and criteria (C) is best matched to 
the methodology of the spectroscopic slit analysis of CIO. 

Close pair fraction predictions for the full uncorrected 
observable Nc and the spectroscopic slit may be found in 
Figure [S] Note that these values are not corrected for line of 
sight projection effects. Here we track objects in the catalog 
from ^ z ^ 3 that have a number density matched to a 
z = 3 sample with Vmax greater than four different critical 
values. By holding a constant Ug cut we examine the changes 
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in Nc at a fixed population size with redshift. The solid black 
curve represents the close pair fraction for a galaxy sample 
with Vin > 133 km s~^ utilising our standard criteria (A). 
The green dodecagonal point and the magenta dash point 
(offset to 2 = 2.9 and z = 3.1 for clarity) are the observed 
serendipitous close pair values of the highest confidence and 
total z — 3 sample, respectively. The solid black triangle 
denotes the results from our fiducial slit criteria (C) and is 
consistent with both observational samples. 

The blue dashed line and the red dot-dashed line repre- 
sent more stringent mass cuts with values of Vi„ ^ 200 km 
s~^ and Vi„ ^ 300 km s~^, respectively, using the standard 
cylinder. These values are included to make a comparison of 
Nc with Ug, the (comoving) number density of galaxies, as 
has been done in work by other authors, even though their 
correlation functions do not closely match the observations. 
In these cases we may examine the evolution of higher mass 
subsamples with a fixed Vin, with redshift. The solid blue 
square and the solid red hexagon represent the expected slit 
model Ncs for these higher mass cuts. As we can see, they 
are both well below the expected Ncs from the LBG obser- 
vations. 

The light blue short-long dashed line and the light blue 
triangle represent the value for our V„ow 5S 142 sample. 
This model also underpredicts Nc, and does so even when 
matched t o the number d ensity of galaxies, Ug, at all red- 
shifts (see iBerrier et al] l|200g ) for more detail). As this 
model demonstrates the same flaws as our Vin model and 
does not reproduce the observed Nc down to low-z, it will 
be excluded from further discussion. 

At our fiducial slit radius of 30.0 ft"^kpc the sample 
impurity is / = 0.1716±0.0111 for the slit and / = 0.2225± 
0.0099 for the full cylinder count. These values imply that 
we have (1 - I) Nc, or a real Ncs ~ 0.042 ± 0.0015 (cf. 
Ncs = 0.0506 ± 0.0017, uncorrected) in the slits and Nc ~ 
0.177 ± 0.005 (cf. iVc = 0.228 ± 0.006, uncorrected) for the 
full cylinder. 

We can generalize the results of Figure [H] to Figure [S] 
We examine the values of Nc in both the spectroscopic slit 
(blue dashed curve) and the full cylinder with radius 10 ^ 
r ^ 30 h~^ kpc (black solid curve) as a function of the 
number density of "galaxies" in the sample, Ug. Here, the 
observational spectroscopic slit results seems too high for 
the observed LBG number density {ug ~ 0.004 ± 0.002 h^ 
Mpc~^) and appears to be more consistent with pair fraction 
at the density of our full Vin ^ 133 km s~^ sample which 
is ~ 5x higher. What causes this inconsistency? We discuss 
this issue in § 13.41 

3.2 Angular Close Pair Fraction 

We also use our simulation to estimate the observed angular, 
or "photometric", close pair fraction, Np. In this case, we 
calculate a pair fraction without using a velocity cut and 
exclude all objects inside the minimum projected radius of 
10h~^ kpc. As discussed above, it is possible to approximate 
the Np at z; = 3 in the simulation box. The simulation is 
120/i~^Mpc long in the comoving coordinates, and thus at 
z — 3 has an approximate redshift length of Sz = 0.18 from 
front to back. Due to clustering effects being less pronounced 
beyond 10/i^^Mpc we may use this box for this purpose, 
despite the redshift range for the survey being Sz ~ 0.6. 
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Figure 8. Expected spectroscopic close pair fractions. Plotted 
are the expectations from 2 = — 3 for different number den- 
sities of galaxies corresponding to the specified velocity cuts at 
2 = 3. The results report the pair fraction in cylinders with radii 
r < 30 h~^kpc (curves). The expected Ncs at 2 = 3 for the 
serendipitous pairs in "observed" spectroscopic slits are also plot- 
ted (points). The solid black curve and triangle point represent 
the close pair fraction of galaxies, Nc, in our fiducial sample of 
haloes using the Vin ^ 133 km s~^ cut at 2 = 3, and the same 
number density at all previous redshifts, matched to the LBG 
correlation function observations. Similarly, the light blue short- 
dash long-dashed curve and empty triangular point represent our 
fiducial Vnow ^ 142 km s~^ sample at 2 = 3. More massive and 
brighter samples are shown with number densities matching the 
values of Vin ^ 200 km s~^ (blue dashed curve, square point) 
and Vin ^ 300 km s^^ (red dot-dashed curve, hexagonal point) 
at 2 = 3. Uncertainties shown for the randomly oriented slits are 
calculated from a combination of jack knife errors and variance 
between the four realisations used. Please note that these points 
are not corrected for sample impurity due to projection effects. 
The green dodecagonal point at 2 = 2.9 and the magenta dash 
point at 2 = 3.1 represent the observed A^c from the high sig- 
nal to noise sample and the full observed sample of 2 ~ 3 LBGs 
l|Cooke et al.ll2010h . respectively, and are offset from 2 = 3 for 
clarity. 



Thus, while we can easily determine Np, we must be cautious 
that this is only an approximate value. 

FigureOprovides us with some estimate of the effects of 
a larger velocity cylinder on the purity of the sample. In the 
case of our photometric sample we are using an equivalent 
velocity difference based on the total length of the box in 
local Hubble fiow as well as the peculiar velocities along 
the line of sight for the galaxies. Similarly we can see from 
Figure |S] that the sample's impurity increases with larger 
sample radius. In the case of 30 ft~^kpc this is ~ 0.20 ±0.01 
from the on the sky dimension of the cylinder alone. 

We examine this solely at 2 = 3 in a fashion similar 
to Figure |9] for our photometric pairs sample. The result is 
Figure [TOl which presents Np as a function of Ug. The solid 
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Figure 9. Variation of the ^ = 3 close companion count (Nc) 
with galaxy number density (rig). The black solid line gives the 
Nc for the full cylinder while the blue dashed line gives the Ncs 
for the randomly oriented spectroscopic slits. The green octagonal 
point represents the observed serendipitous close pair fraction for 
the high confidence sample in the spec troscopic slits, at th e ob- 
served number density, in the sample of lCooke et al.l | |2010| '|. The 
magenta hash is the pair fraction for the total observed sample. 
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Figure 10. The photometric close pair fraction at ^ = 3 as a 
function of the number density of galaxies in the sample. The 
solid black line is the observed photometric close pair fraction. 
The blue dashed line is the close pair fraction for the physically 
associated sample only. This line illustrates only the fraction of 
galaxies that would be mutually inside the same hosting dark 
matter halo. The green octagonal point represents the observed 
close pair fraction for a s ample of photometrically observed LBGs 
from lCooke et al] 1I2OIOI ) (see text for further details). 



black line is the observed A'p, and the blue dashed line is the 
fraction of galaxies that are actually physically associated. 

We find a photometric close pair fraction, N^, of 0.336± 
0.009 before correcting for impurity. After this correction we 
find ~ 0.189, which agrees with our predicted Nc- The 
discrepancy between the "observed" and the physically asso- 
ciated fraction, the "pure" portion of the sample without line 
of sight contaminants, is reasonable due to the large effec- 
tive length of the cylinder used in these measurements. The 
corrected fraction reflects the "real" close (< 30 /i~^kpc) 
companion count and therefore the real fraction of interact- 
ing galaxies. 

These results are interes ting when c o mpare d 
with the results p resent ed in IConseUce et al.l l|2003l ). 



iBertone fc Conselicd (|2009l '). and iBluck et all l|2009l ) which 
produce similar statistics for these high redshift objects. 
The close pair fractions we find are consistent with the 
merger fractions estimated in these works for other ob- 
served large samples of galaxies using close galaxy pairs as 
well as estimates based o n the concentra t ion-a symmetry- 
clumpiness (CAS) method. IConselice et al.l l|2003h estimates 
an apparent merger fractions of bright LBG's at z ^ 2.5 to 
be between 40 — 50 per cent. iBluck et al.l (|2009| ) identifies 
82 massive galaxies, M« > lO^^'^M© in a redshift range 
of 1.7 — 3.0. These are further divided into a sample of 
44 galaxies between 1.7 < z < 2.3 and 38 galaxies from 
2.3 < z < 3.0. Close pair fractions are estimated for the 
two samples by identifying all imaged galaxies within 
±1.5 magnitudes of the sample galaxies magnitude that 
reside within an Rmax < 30 kpc (physical, h = 0.7) and 



statistically corrected for impurities. Please note that 
only the original 82 galaxies have redshift information, 
the rest of the galaxies included in the measurement are 
photometric pairs. The observed Np for the 1.7 < z < 2.3 
(2.3 < z < 3.0) are Np = 0.19 ± 0.07 {Np = 0.40 ± 0.10). 
These observed high redshift values are consistent with our 
estimates in both th e corr ected and uncorrected sample. 
iBertone fc Conselicd (|2009l ) examines models of galaxy 
merger rates comparing simulations and observational 
results. This work also finds a high merger rate at high 
redshift. 



3.3 Comparison To Previous Results 

In order to demonstrate the usefulness of this technique in 
calculating Nc across a wide range of redshifts we make a 
comparison to several previous existing measurements at low 
redshift (see Berrier et al. (2006) for further discussion on 
these samples). The samples used in this comparison are ex- 
tracted from t he SSRS2, CNOC2, an d the DEEP2 surveys 
iPatton et al] (|2002l ): iLin et all (|2004l ). These surveys pro- 
vide several candidate definitions for a close pair. To match 
the technique used in the observations we use criteria (A). 
In order to determine the sample of haloes used in our sim- 
ulation, both host dark matter haloes and substructure, we 
match the number density of galaxies from the observations 
to the number d ensity of haloes in t he simulation using the 
Vin model as in iBerrier et al.l (|2006l ). 

To extend this method to 2; ~ 3, we examine observed 
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colour-selected LBG pairs with 10 — 30h~^ kpc separations 
from our LBG survey and the larger survey of SOS and find 
an observed photometric pair fraction of Np = 0.047±0.035. 
This result is corrected for spatial impurities that are esti- 
mated in a manner similar to that for the low-z samples. 
We estimate the impurity using mock catalogs constructed 
to the exact field dimensions and number densities of each 
observed field. We then distribute mock galaxies using red- 
shift distributions and interloper fractions determined by the 
photometric selection function. The number of close pairs 
observed in projection is then corrected to align with the 
fraction that is found to consist of true pairs in three di- 
mensions. 

We then calculate the z = 3 data in the simulation in 
the same manner as the low-redshift data and matched it to 
the LBG number density. From our definition of the close 
pair fraction, a decrease in number density corresponds to a 
similar decrease in the close pair fraction. If LBGs randomly 
comprise ~ 1/5 the number of massive haloes as the number 
densities imply, and thus approximately 1/25 of our sample 
would be composed of LBG-LBG pairs, we would expect a 
lower limit of Nc ~ 0.228 x (1/5), or ~ 0.0456 {Nc ~ 0.177 x 
(1/5), or ^ 0.0373 for the purity corrected sample) for the 
observed LBG-LBG pair fraction. Our simulated value is 
stiU larger than this at Nc = 0.083 ±0.005 {Nc ~ 0.065 with 
purity corrections). For comparison, we look at Np for our 
photometric pair value after corrections for impurity, as this 
is the only estimate of the full Nc available. This assumes 
that all galaxies in the fiducial simulation sample which are 
within the same host halo will be observed, irregardless of 
the rig issues discussed in § 12.81 

The results of this comparison are presented in Fig- 
ure 111! Here the solid black line, black points, and solid 
square points represent the pair fractions calculated from 
our Ug matched samples from z = to 2 ~ 3. The empty 
squares are observe d values from the DEEP2 Survey taken 
from fields 1 and 4, iLin et al] (120041 ) and the triangles are 
from the CNOC2 and SSRS2 survevs. IPatton e t al. (2002). 
The value of Nc from the z 3 observations is indicated by 
the green dodecagonal point at 2; = 3. The full fiducial sam- 
ple Nc for 2 = 3 is shown connected by the dashed line. The 
Nc values determined from the simulations are closely com- 
parable to the matched-density observations in nearby red- 
shift bins across all redshifts observe d. The technique us ed 
here is more completely described in iBerrier" et all (|2006l ). 

3.4 Galaxy Number Density Problem Revisited 

We have matched the 3-D two-point correlation function 
to produce our fiducial sample of high redshift galaxies. We 
find that if LBGs comprise all massive haloes in the matched 
simulation sample, this produces a number density that is 
too high by a factor of ~ 4.75 when compared with the 
observed LBG number density, as has been similarly found 
by other authors. In addition, this sample produces a value 
of Np and Nc in the cylinders that is too high by a similar 
factor but a value of Ncs in the spectroscopic slits (that are 
biased to detecting pairs with separations of <~ 20 h~^ 
kpc because of their geometry) that is consistent within the 
errors of the observations. 

In contrast, forcing the sample to match the observed 
LBG number density yields an Np and cylinder Nc similar 
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Figure 11. Comparison of Nc from 2: = to 2: = 3 in the 
simulation with measurements from surveys. The black points, 
solid squares, and solid line represent the Nc calculated from our 
Ug matched samples from z = to z = 3. The dashed black line 
represents Nc for our fiducial sample at 2: = 3. The empty blue 
squares a re data points f rom the DEEP2 Survey taken from fields 
1 and 4, l|Lin et al.ll2q04h. The hollow re d triangles are data from 
CNOC2 and SSRS2. jPatton et al.ll2002l) . The green dodecagonal 
point at 2 = 3 is the Np measured from the LBG surveys of 
.Cooke et al.l (l2005h and lSteidel et al.l (120031 '). 



to the observations but an Ncs that is significantly lower 
than that observed in the slits. Moreover, such a sample 
{Vin > 200 kms~^) results in a poor fit to the observed LBG 
correlation function and corresponds to haloes too massive 
to be reconciled with the mass of LBGs from clustering anal- 
ysis. 

LBGs do not comprise all massive galaxies at 2 ~ 3. 
Other identified populations include sub-mm galaxies, pas- 
sive and star forming BzK galaxies, DRGs, and other galaxy 
types with typically lower UV luminosities then LBGs. If 
LBGs represent ~ 1/4.75 of the matched massive haloes 
in our 2 ~ 3 sample, then the agreement in number den- 
sity produces an A'j, and cylinder Nc that are in very good 
agreement with the observations but an Ncs in the slits that 
is significantly lower than the observations. In this case, the 
sample still matches the LBG correlation function (as they 
are pulled from the same parent population) and thus cor- 
responds to haloes with the same mass as the observations. 

Dust obscured galaxies may make up a fraction of the 
"missing" massive haloes, however, infrared and sub-mm 
surveys recover only a small fraction of the number necessary 
to reconcile the difference. The latter interpretation more 
accurately refiects the observed LBG population and leaves 
only the < 20h~^ kpc physical pair fraction in disagreement. 
The small-scale behaviour (one-halo term regime) measured 
in accurate 2-D correlation functions that utilise deep, wide- 
field imaging helps provide the solution to the remaining 
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disagreement and offers interesting insight into the nature 
and detectability of LBGs. 



3.5 Small-scale behaviour: The one-halo term 

Measurements of the angular (2-D) correlation function 
(ACF) at z ~ 4 in the Subaru /X MM-Newton Deep Field 
over 1 deg-^ (|Ouchi et all |2005| ) and at z ~ 3 in the 
Canada-France-Haw aii Telescope Legac y Survey 4 square- 
degree Deep fields (|Cooke et alj 12012 1) are able to utilise 
a large number (^10"* — 10^) of LBGs to accurately probe 
the ACF down to mu ^ 27 from relatively small to large 
scales (~ 0.5 — 10 Mpc, comoving). Both efforts witness 
a distinct break in the form of the ACF at small scales from 
the power law fit over larger scales that may provide insight 
into the discrepancy in the observed and expected close pair 
fractions in the spectroscopic slits. 

In order to measure the angular correlation function at 
z = 3 in the simulatio n, we will follow the method used in 
previous works such as lConrov et al.l (|2006| ) and utilize the 
Limber transformation. 



10 



(4) 



where r is the comoving distance at z and N{z) is the nor- 
malized redshift distribution of the galaxies in the observed 
sample. 

Figure [12] presents the ACFs of our simulation sam- 
ples and the two observational datasets. As demonstrated 
in Figure [121 the form of the fiducial {Vi„ > 133 km s~^) 
and observational ACFs in the outer regions are consistent. 
However, the inner regions show a marked discrepancy. Re- 
gardless of the number density, the standard technique of 
halo assignment cannot reproduce the features of the LBG 
angular correlation function on both large and small scales. 
In addition, the mismatch cannot be corrected by any scaling 
of the data via an integral constraint. Our standard abun- 
dance matching model is unable to reproduce the observed 
break from power law behaviour near 150 - 200 h-i kpc m 
comoving coordinates to match both the one-halo and two- 
halo components of the correlation function to the accuracy 
of the data without incorporating assumptions of the form 
of the Limber equation in the inner region. 

If we follow a standard subhalo abundance matching 
scheme we observe a Nc in the cylinder and A'p from our 
mock photometric sample which are consistent with the ob- 
served pair counts, however, we do not produce the correct 
Ncs for the mock spectroscopic slits (pairs at < 80 /i^^kpc 
in comoving coordinates, or < 20 /i^^kpc in physical coordi- 
nates). Nc is smaller than the observed fraction by a factor 
consistent with the decrease in amplitude of the one-halo 
term in the correlation function over the same separations. 

To furt her test this result, we analy se the Millennium 2 
simulation ( Bovla n-Kolchin et al.l [20091 ). The Millennium-2 
simulation has five times the spatial resolution of the origi- 
nal Millennium simulation, in this case a Plummer equiva- 
lent softening of 1 /i~^kpc. This resolution is more than ade- 
quate for our needs. The sample we select from Millennium-2 
is well above the resolution limits of the simulation, and thus 
provide a further test that our results are not adversely ef- 
fected by resolution issues. Using this simulation, we find the 
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Figure 12. The angular correlation function for samples of dark 
matter haloes pulled from our primary sample. Two subsamples 
are taken from the dark matter halo catalog using different Vin 
cuts and are compared with two observational sets. The solid 
black line is our fiducial sam ple which matches the t hree dimen- 
sional correlation function of lAdelberger et al.l ([200^ . The green 
long-dashed line is our best match to data from the Millennium-2 
simulation sample with Vin ^ 123 km s~^. The magenta dot-short 
dashed line and cyan short-long dashed line represent abundance 
matched samples with Vin ^ 180 km s~^ from Millennium-2 and 
Vin ^ 200 km s~^ from the Z05 simulations respectively. The 
red short dashed sample repr esents the observed a ngular correla- 
tion function for 2 ~ 3 from ICooke et al. I 1I2OI2 I) , and the blue 
dot-long dashed lin e illustrates a sample of LBGs at 2 = 4 from 
lOuchi et al.l ||2005| ). Both observational samples demonstrate a 
strong break from a single power law around 100 h~^ kpc. None of 
the simulation samples are able to reproduce this feature through 
our standard abundance matching techniques. The two dark red 
solid vertical lines represent the range, in comoving coordinates, 
that we identify galaxy- galaxy pairs. 



closest match to the two point correlation function of the ob- 
servations to be a cut of Vin ^ 123 km s~^. Haloes with this 
criterion show a good agreement to the data, with a power 
law fit of ro = 3.98 ± 1.44 Mpc h'^ and 7 = -1.47 ± 0.27. 
In addition, this sample has a number density rig = 0.017 
h'^ Mpc^'^, approximately 4.25 times larger than the ob- 
served of LBGs and similar to the overdensity of our pri- 
mary simulation sample. We have included the Millennium-2 
ACFs for the abundance matched sample and for the 3-D 
correlation function matched sample in Figure 1121 Again, 
the simple abundance matching techniques are not capable 
of reproducing the shape of the angular correlation function. 

This selected sample produces a, Nc — 0.1521 ± 0.0006 
for our fiducial criteria (A) in the cylinder, and Ncs = 
0.033 ± 0.001 for our fiducial criteria (C) in the slits. The 
sample shows an impurity of 0.1433 ± 0.0176 in the cylinder 
and 0.1167 ± 0.0202 for the slit, thus our corrected values 
are 0.1303 ±0.0059 and 0.02919 ±0.0014 for the fuU cylinder 
and the spectroscopic slit respectively. For our line of sight 



© 2012 RAS, MNRAS OOO.fTHTSl 



14 Berrier & Cooke 



mock photometric sample we find A'p = 0.2299 ± 0.0144, 
before corrections for impurity and 0.1370 ±0.0154 after. 

To match the number density of observed LBGs we 
would select a sample of haloes with Vi„ ^ 180 km s~^. As 
with the previous number density matched sample this has a 
power law fit of ro = 5.17 ±2.24 and 7 = -1.56 ±0.30. This 
sample produces a close companion count of only Nca — 
0.0101 ± 0.0023 in the spectroscopic slit, Ncs = 0.0089 ± 
0.0022 after impurity corrections, iV^ = 0.0505 ± 0.0060 for 
the full cylinder, only 0.0442 ± 0.0061 after being corrected 
for impurity, and A'p = 0.0682 ± 0.0092 photometrically, 
0.0471 ± 0.0075 after correcting for sample impurity. 

All results for our primary samples in the Z05 simula- 
tions and in Millennium-2 from here and from § 13.11 § 13.21 
are summarised in Table [T] 

3.6 Discussion 

We find that matching haloes in our simulation to the ob- 
served LBG number density or the LBG 3-D correlation 
function and mass using a simple prescription can generate 
informative close pair statistics. We find that the low density 
simulation samples are able to reproduce the total observed 
Nc and Np in cylinders, but underpredict the fraction ob- 
served serendipitously in spectroscopic slits. In contrast, the 
higher density 3-D correlation function matched sample is 
able to reproduce the spectroscopic slit fraction, but over- 
predicts the number of observed galaxy-galaxy pairs. Our 
numerical/analytical simulation is not resolution limited in 
this sense and the discrepancy at small scales occurs above 
the resolution limit of the Millennium-2 simulation. At these 
separations, many of the luminous galaxies sharing these 
haloes are either interactions or imminent interactions. This 
finding provides an interesting avenue to quantify the spatial 
behaviour of LBGs and sub-halo assignment schemes. 

Galaxy interactions can generate a significant enhance- 
ment in their luminosit ie s from the close int e ractions (e.g. 
Larson fc TinslevI Il978l: iBarton et al.l I2OO0I : lEllison et at] 
20081 : [Bridge et al.ll201oi : IWong et al.ll201lD . Moreover, it is 
possible that the luminosity enhancement at high redshift is 
equivalent to, or higher than, that observed at low redshift 
as a result of the higher gas fractions in LBGs. 

The LyQ emission versus separation relationship of ob- 
served close LBG pairs found in CIO supports this picture. 
All of the spectroscopic < 20 /i~^kpc physical close pairs, 
and thus interacting systems, exhibit Lya emission as com- 
pared to ~ 50 per cent of the full population. A fraction 
of the observed Lya emission of each interacting galaxy is 
likely to be a signature of enhanced star formation. This be- 
haviour may extend to the Lya emitter (LAE) population 
as well (see CIO). 

Lower-luminosity LBGs t ypically have lower m asses 
(|Giavalisco fc DickinsonI I2OOII : iKashikawa et al.l I2OO6I ) and 
the higher density of these haloes yields a higher interaction 
fraction as compared to our fiducial (m_R < 25.5) sample. 
The luminosity enhancement from interactions would boost 
a fraction of lower-luminosity LBGs above the magnitude 
selection cut-off of our sample. This process would create an 
increase in the number of < 20 h~^kpc physical close pairs 
detected in the observations that are not represented in the 
simulation analysis. 

Our adopted halo assignment (§2.3), does not account 



for enhanced star formation. Moreover our Vnow sample 
(§ 12. 3[) models haloes where baryons are stripped during 
infall. Such galaxies would have a decrease in luminosity as 
compared to our Vi„ model, and we find the V„ow model 
predicts fewer close pairs detected in the slits. This result 
implies that a model which instead includes an appropriate 
luminosity enhancement per baryon for infalling haloes over 
the standard assumptions of abundance matching will pre- 
dict a higher fraction of close pairs in the slits as is seen in 
the observations. 

Our results and the proposed scenario remain con- 
sistent within high redshift measurements and frac- 
tions of close or interacting/merging LBGs by vari- 
ous means when considering the samples studied (c.g^ 
Conselice ct al. 2003; Lotz ct al. 2006; Bcrtonc & Consclica 
20091: (F5rster Schreiber et al.l I2OO9I : iBluck et all 1200a 
Law et al.ll20ldV_ In addition , obse r vations of low redshift 



LBG analogs (|Overzier et all |2009| . I2OI0I : iGoncalves et all 
l2010h . which are matched to LBGs in essentially every way 
(stellar mass, gas fraction, star formation rate, metallicity, 
dust extinction, physical size, gas velocity dispersion, etc.), 
show from optical imaging that the bulk of these objects 
are undergoing interactions even though the UV imaging is 
inconclusive. Simulated to high redshift, these objects are 
consistent with the properties and observations of 2: ~ 3 
LBGs^ 

iLaw et"ai] l|2012l ) uses Hubble Space Telescope imaging 
in the restframe optical to estimate the number of real close 
pairs in a sample of galaxies observed between 2.5 ^ 2 ^ 3.6. 
Galaxies with spectroscopically determined redshifts and 
magnitudes between H — 22.0 — 24.0 were compared to ob- 
jects within 5 — 16 h~^ kpc with no more then one magnitude 
difference. These results were statistically corrected for false 
close pairs and produce a value of Nc = 0.17^0 ^8 ^'^^ 2: ~ 3 
LBGs. This close pair fraction is consistent with our results. 

In our fiducial model, selected by matching the two 
point correlation function, we have not truly required all 
galaxies to be visible either due to dust extinction or low 
star formation, thus we do not have to match the observed 
number density a s we would in Sub H alo Abundance Match- 
ing. The work of iReddv et al.l l|2008l ) suggests that the rest 
frame UV luminosity of galaxies at these redshifts are typi- 
cally extincted by a factor of 4 — 5 in fiux. Our typical halo 
masses at accretion are M ^ 10^^'^*. All but ~ 1 per cent 
of our halos have a mass at accretion above the minimum 
mass oi Mmin ^ xg^^'^*'''^ h~^ Mq su ggested by ha l o mod - 
eling for z ~ 2 star forming galaxies in lConrov et al.l l|2008h . 
where the number density of halos does not match the ob- 
served rig for that population. 

Our number-density matched sample implies that every 
halo hosts a luminous LBG. In this case, the lower fraction 
of close pairs at < 20 kpc separations as compared to the ob- 
servations goes against interacting galaxy behaviour. More- 
over, we know from high redshift surveys that LBGs com- 
prise a large fraction, but not all, detectable galaxies at high 
redshift and that haloes indeed host massive galaxies that 
are not detectable using Lyman break techniques. Thus, we 
are forced to consider higher densities samples. In order to 
generate 4.75 times the observed LBG density of the corre- 
lation function matched sample, we need to integrate down 
the faint-end of t he luminosity function below R = 25.5 
by ~1 magnitude jSawicki fc Thompsonll200"^ : iReddv et al.l 
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Table 1. Nc and Np for all samples 



kJllllLilOiLlLJll Oo i 1 Lt-i lex 


N 






\^\-JL J. CL^ cs 


rig 


Z05 (A) 


0.228 ± O.OUd 


0.177 ± 0.005 


0.0293 ± 0.0013 


0.0233 ± 0.0011 


0.019 


Z05 (B) 


0.246 ± 0.006 


0.193 ±0.005 


0.0375 ± 0.0016 


0.0305 ±0.0014 


0.019 


Z05 (C) 


0.275 ±0.007 


0.218 ±0.006 


0.0506 ±0.0017 


0.0419 ±0.0015 


0.019 


Z05(R)(A) 


N/A 


N/A 


0.0220 ±0.0234 


N/A 


0.019 


Z05(R)(B) 


N/A 


N/A 


0.0300 ± 0.0178 


N/A 


0.019 


Z05(R)(C) 


N/A 


N/A 


0.0430 ± 0.0145 


N/A 


0.019 


Z05(Pliotometric) 


0.336 ± 0.009 


0.189 ±0.003 


N/A 


N/A 


0.019 


Millennium-2 (A/C) 


0.1521 ±0.0006 


0.1303 ±0.0059 


0.033 ±0.001 


0.02919 ±0.0014 


0.017 


Millennium-2 (Photometric) 


0.2299 ±0.0144 


0.1370 ±0.0154 


N/A 


N/A 


0.017 


Observations Photometric 


0.047 ± 0.035 


N/A 


N/A 


N/A 


0.004 


Observations Spectroscopic 


N/A 


N/A 


0.047 ±0.015 


N/A 


0.004 



The Z05 models are our standard simulations. Millennium-2 simulation values are reported using criteria (A) in the cylinder and 
criteria (C) in the slit. The number density of the abundance matched sample used in Figures 8 and 11 are rtg = 0.004 h^ Mpc-3 



1200^ '). Although a fraction of the higher luminosity haloes 
in this magnitude range are expected to have sufficient star 
formation enhancement to enter into our magnitude cut, the 
lack of knowledge of the typical enhancement in the FUV 
at z ~ 3 makes it unclear if such a fraction is large enough 
to match that observed in the spectroscopic shts. If a signif- 
icant fraction of the star forming galax ies are obscured by 
dust as suggested bv lLacev et al.l (|201lh . a random subsam- 
ple would have approximately the same correlation function 
and produce similar close pair fractions. As a result, the 
enhanced fraction of lower-luminosity galaxies may include 
high star formation rate, dusty galaxies that may experience 
a larger magnitude increase from the effects of interaction 
and morphological disruption. 

Finally, assuming that LBGs co mprise >50% of all 
star forming galaxies at z ~ 3 (|Reddv et al.l l2005l : 
iMarchesini et all 120071 ') . the results of the correlation func- 
tion matched sample suggest that after considering dust- 
obscured galaxies, ~3/5 of haloes are not observed us- 
ing any high redshift detection technique. Driven by (1) 
the power of the sim ulation to predict t he close pair frac- 
tion down to z = (|Berrier et al.l [20061 ). (2) the evidence 
that a bundance m atching may be used with high redshift 
LBGs ijConrov et al. 2006). (3) the possible direct correla- 
tion between halo mass and UV luminosity at this epoch 
ijSimha et al.ll2012l : IConrov fc Wechsieill2009l ). and (4) the 
equivalent overdensity of similarly matched massive haloes 
in other simula tions including si mulations using different 
approaches (e.g lLacev et alJl201ll l. the unaccounted z = 3 
haloes in the correlation function matched sample likely re- 
flect a similar number of real haloes in the Universe. If true, 
these haloes must either have highly obscured star formation 
that is not detected by curren t high redshift selection tech- 
niques (e.g. iLacev et al.ll201 such as IR and sub-mm sur- 
veys, or they must be massive galaxies with inherently low 
star formation rates that are below the detection thresholds 
of current facilities. It may be the case that galaxy interac- 
tions is the cause for the initial starburst, or "turn-on", of 
many of these undetected haloes within our sample. Thus, 
a combination of all of the above affects resulting from in- 
teraction induced star formation may provide a plausible 
explanation for the larger fraction of observed (< 20 kpc) 
serendipitous pairs in spectroscopic slits when assuming 



LBGs comprise ~l/5 of the correlation function and mass 
matched sample and ~l/5 the reported close pair fractions. 
One means to probe such massive haloes independently of 
their luminosity is via quasar absorption line systems , in par- 
ticular, the ubiquitous damped Lya systems (.Wolfe et al.l 
2005), which have been shown to be as sociated with massi ve 
systems that cluster s i milar to LBGs jCooke et akl l2006al lbl: 
iNagamine et al" I I2OO6I : iLee et al.|[2oTll ). We are engaged in 
investigation that is testing various components of this sce- 
nario on several fronts. 



4 CONCLUSIONS 

We have matched the 3-D two point correlation function of 
a sample of 2 ~ 3 LBGs to a sample of haloes from our pri- 
mary numerical/analytic cosmological simulation. Using this 
sample we have mocked observations of simulated spectro- 
scopic slits and of photometric observations of these galax- 
ies. We also test our model with data from the Millennium-2 
simulation for verification of our results. This work has led 
us to several interesting results which we summarise in the 
points below, see also Table [T] above. 

• We demonstrate that neither standard SubHalo Abun- 
dance Matching (SHAM) or a two-point correlation function 
and mass matching scheme completely reproduces the obser- 
vational results. Neither model can reproduce galaxy clus- 
tering features and Ug at same time. Explicitly, we find that 
the standard SHAM does not reproduce the serendipitously 
observed Ncs , and the break in the LBG correlation function 
at very small scales (< 20 h~^kpc physical, ~ 80 h~^kpc co- 
moving) . 

• The number density of our candidate LBG sample is 
~4.75 times the observed LBG number density. The im- 
plication is that only i halos above a fixed mass are de- 
tectable LBGs. These re s ults a r e con si stent with the re - 
sults of iNagamine et all tOOA [200^ ; iLee et al.l (|201lD : 
iDave et al. l (|2000l ): lLacev et al.ll|201lD and others which find 
a similar overdensity using other types of simulations. 

• We find an observed close pair fraction Nc = 0.228 ± 
0.006, which implies an impurity corrected close pair fraction 
of iVe = 0.177 ± 0.005, (~ 18 per c ent) This result is con - 
sistent with the previous results of iGonselice et all (|2003h : 
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iBertone" 



I 2012h 



& Conselice ('2009'); 'Bluck et al.' ('2009 1'): iLaw et alj 
Lotz et al (.2006 ); Fors tcr Schrcibcr et al.l (l2009f ) 



considering the uncertainties specific to those studies and 
pair fraction/merger rate assumptions. 

• Our simulated matched spectroscopic shts produce a 
close pair fraction of Ncs ~ 0.0506 ± 0.0017 for our fiducial 
sample, defined to have a maximum velocity separation of 
±500 km s~^and an on the sky separation of ^ 30.0 h'^kpc. 
This is similar to the observed fraction of serendipitous spec- 
troscopic close pairs of Ncs = 0.047 ± 0.015 for the full ob- 
served LBG sample and Ncs ~ 0.071 ± 0.023 for the highest 
signal to noise ratio sample. After correcting for false close 
pairs which may be observed we find Nc = 0.0419 ± 0.0015 
in the simulation. 

• If we examine our catalogs using randomly selected 
slitlets to generate a more generic result, we find that the 
expected fraction of LBG pairs that fall serendipitously into 
the slitlets to be Ncs = 0.0430 ± 0.015 before correcting for 
sample impurity. 

• We find a photometric close pair fraction of Np — 
0.336 ± 0.009 and after correcting for the sample impurity 
we find A'p ~ 0.189. The latter fraction reflects the "real" 
number of close pairs and therefore the real number of po- 
tentially interacting galaxy pairs. As mentioned above, only 
a portion of the corrected value will be observable LBGs. 
The difference in Ug is a factor of ~ 4.75 leading to a cor- 
rected value of A'p ~ 3.99 per cent, which is consistent with 
our photometric pair fractions Np = 0.047^:0.035 estim ated 
from our survey and the survey of lSteidel etahl (|2003| ). 

• The analysis of the sample taken from Millennium- 
2 produces similar results to our primary simulation. In 
Millennium- 2, we find the correlation function-matched 
sample to be overdense in comparison with the observa- 
tions by a factor of 4.25. This selected sample produces a 
Nc = 0.1521±0.0006 in the cylinder, and Ncs = 0.033±0.001 
in the slits. The sample shows an impurity of 7 = 0.1433 ± 
0.0176 in the cylinder and I = 0.1167 ± 0.0202 for the slit, 
producing corrected values of Nc — 0.1303 ± 0.0059 and 
Ncs = 0.02919 ±0.0014 for the fuU cylinder and the spectro- 
scopic slit respectively. For our line of sight mock photomet- 
ric sample we find Np = 0.2299±0.0144, without corrections 
for impurity and A'p = 0.1370 ± 0.0154 after correction. 

The excess of close (interacting) pairs 20 h"'^ kpc 
(physical) and the inability for the standard abundance 
matching with monotonic UV-mass halo assignment to de- 
scribe the steep slope in the observed LBG correlation func- 
tion at very small scales provides insight into triggered star 
formation and the detectability of LBGs (and LAEs) at 
z = 3. Our results imply that the spectroscopic slit close 
pair fraction and the break in the correlation function repre- 
sent the detection of either a fraction of less massive (higher 
density) LBGs with luminosities below our magnitude cut 
(m_R = 25.5) as a result of an enhancement in luminosity 
from interactions, the "turn-on" of massive haloes with pre- 
vious low star formation as a result of interaction, or, likely, 
a combination of both cases. 

We find that LBGs likely represent ~ 20 — 25 per cent of 
all massive {Vi„ > 133 km sec~^) haloes at 2 ~ 3 based on 
the results of the analysis of our simulation, the Millennium- 
2 simulation, simulation analyses by several other authors, 
and the power of our simulation analysis to predict the close 



pair fraction from 2 = 1 to 2 = 0. The full census of detected 
star forming galaxies selected by various criteria suggest that 
LBGs likely account for > 50 per cent of the massive haloes 
at 2 ~ 3. The remaining fraction is likely populated by sys- 
tems with low star formation rates and/or systems that are 
not detected using current selection techniques. DLAs are a 
promising means to explore the remaining fraction of mas- 
sive haloes because they probe galaxy haloes randomly, in- 
dependent of luminosity, they have a high nu mber density, 



and they are found t o resid e in massive haloes (ICooke et al 
2006b; F vnbo et al .' 2003, 2008, 201(^, I2OIII : iMoUer et al 
2002 : MoUer et al.ll2004 : ,Schave,2001. ) . 



The statistics generated from our mock spectroscopic 
slits with the serendipitously confirmed close pairs from ob- 
servations provides a potentially powerful tool to estimate 
the behaviour and nature of LBGs and the enhanced star 
formation rate from LBG interactions. 
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